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NON EQUILIBRIUM STATIONARY STATE FOR THE SEP WITH 

BIRTHS AND DEATHS 
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o 

CN ' Abstract. Wc consider the symmetric simple exclusion process in the interval Ajv := 

^H . [— A^, A^] n Z with births and deaths taking place respectively on suitable boundary intervals 

/+ and /_, as introduced in De Masi et al. (J. Stat. Phys. 2011). We study the stationary 
measure density profile in the limit A^ — > oo. 



cd 



OO 

^'. 1. Introduction 

♦-^. 

r-| ' This paper is a follow-up of the study initiated in [1] , [2] , where current reservoirs in the con- 

^ ■ text of stochastic interacting particle systems have been proposed as a method to investigate 

stationary non-equilibrium states with steady currents produced by action at the boundary. 
Due to the particular difficulties in implementing this new method, we consider the simplest 
possible particle system. The bulk dynamics is the symmetric simple exclusion process 
^ ■ (SSEP) in the interval Ajv = [— A^, A^] fl Z (A^ a positive integer and N ^ oo eventually), 

— ? ! namely the state space is {0, l}^'^ (at most one particle per site): independently each particle 

t^^ I tries to jump at rate A^^/2 to each one of its nearest neighbor (n.n.) sites, the jump then 

takes place if and only if the chosen site is empty, jumps outside A^v are suppressed. To 
^ ■ induce a current we send in particles from the right and take them out from the left, and 

^ ! would like this to happen at rate Nj /2, j > a fixed parameter independent of A^. Due 

to the restrictions imposed by the configurational space, we have to be more precise when 
defining this dynamics. For this we fix a parameter K > 1 (an integer) and two intervals J± 
kj ■ of length K at the boundaries: I+ = [N - K +1,N] and /_ = [-A^, -A^ + K - 1]. At rate 

^ ! Nj /2, when /+ is not totally occupied, we create a particle at its rightmost empty site; with 

the same rate, unless /_ is empty, we take out a particle from its leftmost occupied site. In 
case /+ is already full, or /_ empty, the corresponding mechanism aborts. 

In [1], [2] we have proved that at any time t > propagation of chaos holds and that in the 
limit A^ — 7- OO the hydrodynamical equation is the linear heat equation: 

|P('^'^) = ^|:^P('^'^)' rG(-l,l),t>0, 

p(r,0)=po(r), p(±l,t)=M±(t), (1.1) 

where po(') is given but u±{t) are solutions of a nonlinear system of two integral equations, 
see (2.5) below. 

The goal of this paper is to investigate the limiting density profile (as A^ — )■ oo) of the (unique) 
invariant measure of the process. The main result is Theorem 2.2, which shows that this 
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rescaled limiting profile coincides with the unique stationary solution of (1.1). In particular, 
taking into account the validity of the Fourier law, proven as Theorem 2 in [1] , we see that 
the effective current in the stationary regime is strictly smaller than its desired maximum 
value which is min{j'/2, 1/4}, but this value is indeed approached by letting K — > oo. 



2. Model and main results 

Particle configurations are elements rj of {0, l}^'^, rj{x) = 0, 1 being the occupation number 
at X G A^r. We consider the Markov process on {0, l}^'^ defined via the generator 

where L^ = L^ + + L^, „ and 

x=-N 



(2.1) 
V D^TiixMfin'^''^^- f(n) . 



Lt,±f{v) := ^ E D±v{x)[fiv^''^) - fiv) 



xei± 

7]^^^ being the configuration obtained from 77 by changing the occupation number at x, r^(^'^+i) 
by exchanging the occupation numbers at x, x + 1; for any u : K^ — )■ [0, 1] 

D+u{x) = [1 - u{x)]u{x + l)u{x + 2) . . . u{N), x e 1+ 

D^u{x) = u{x)[l-u{x-l)][l-u{x-2)]...[l-u{-N)], x G /_. (2.2) 

Given po ^ C'([— 1, 1], [0, 1]), let v^^^ be the product probability measure on {0, l}'^^ such 
that iy^^'{ri{x)) = po{N~^x) for all x G Ajy. Let P,^(iv) denote the law of the process with 
initial distribution u^^'^ and Ej,(]v) the corresponding expectation.^ 

The following theorem has been proven (in a stronger form) in [1], [2]. The statement below 
contains all what is needed in the present paper. In the following, for n a positive integer we 
write A^^ for the set of all sequences (xi, ...,Xn) in A^ such that Xj 7^ Xj whenever i ^ j. 

Theorem 2.1. There exists t > so that for any po os above and any n > 1, 



lim 



E^^M){Y[vix^,t))-Y[E^^N,{r]ix,,t)) =0, foranyt<T\ogN. (2.3) 



j=l i=l 



Omitting the initial profile to avoid too heavy notation. 
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Furthermore 

lim sup sup |Ey(jv) (?7(x, t)) — p(A^~'^x, t)| = 0, (2.4) 

where the function p{r,t) solves the heat equation -^ = 2d^' ^ ^ i~^^ l),t > with initial 
datum pQ and boundary conditions p(±l,t) = u±{t), the pair {u^{t),u_{t)) being the unique 
solution of the non linear system 

u±{t)= [ Pt{±l,r)po{r)dr + ^- [ \p,{±l,l) {l - u+{t - sf) 
J[-i,i] ^ Jo "- 

-P,(±l, -1) (1 - (1 - M„(t - s))^) }rfs, (2.5) 

where Pt{r,r') is the density kernel of the semigroup (also denoted as Pt) with generator 
A/2, A the laplacian in [—1, 1] with reflecting, Neumann, boundary conditions. 
The function p{r, t) satisfies 

^^|.=i = j(l - u^{tr\ ^^|.=_i = J(l - (1 - n^(t)r). (2.6) 

or or 

Remark. The following is the integral form of the macroscopic equation: 

p{r,t)= [ p,{r,r')p{r',0)dr' + ^- [ \p,{r,l) {l - p{l,t - sY') 

J[-i,i] ^ Jo ^ 

-Psir, -1) (1 - (1 - p(-l, t - s)f) }ds. (2.7) 

It will be convenient to recall the expression for the density kernel Pt{r,r') in terms of the 
Gaussian kernel 

Gt{r, r') = ^= , r, r' G R, (2.8) 

V 27rt 
as 

Pt{r,r') = Y^ Gtir,r") for rV ±1 
P,(r,±l) = Yl 2G,(r,r"), (2.9) 

r":tp{r")=±l 

where ■?/; : M — )■ [—1, 1] denotes the usual reflection map: ip{x) = x for x G [—1, 1], ip{x) = 
2 — X for X G [1, 3], ^ extended to the whole line as periodic of period 4. 

Notation. Ptg{r) = J Pt{r,r')g{r')dr', for g a bounded continuous function, t > 0. 



The main result of this paper is about the density profile of the unique invariant measure 
Pn- 
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Theorem 2.2. For any integer k > 1 we have 



lim max 



f^^{ri{x,)---ri{xk))-p*{xi/N)---p*{xk/N) =0 (2.10) 



where p*{r) is the unique stationary solution of the macroscopic equation. Namely p*{r) = 
Jr + l 

J = j(l — a ), with a the solution of a{l + ja ~ ) = j -\ — • (2-11) 

By Theorem 2.2 it follows that fiN concentrates on a L^-neighborhood of the limit profile 
p*: let r G (0, 1) and 



Then for any a G (0, 1) 



..X) 

x£Ajf:\x-rN\<£ 



lim 



^N[f_\p^''^\r:ri) - P*{r)\dr) =Q 



Theorem 2.2 will follow from • uniformly on the initial datum po the solution p(r, t|po) of 
the macroscopic equation (2.7) converges in sup norm to p* exponentially fast, see Theorem 
4.1 below; • for any integer /c > 1, 



K K 

E^[l[r]{xi,t)^ -l[p*{x,/N)\ = (2.12 



lim lim max max 

We are also working on an extension of the theorem where we prove exponential convergence 
in time to p^r uniformly in A^. 



3. MONOTONICITY PROPERTIES 

We consider the space {0, l}^'^ endowed with the usual partial order, namely we say that 
?7 < ^ iff r]{x) < ^{x) for all x G A^r. The following proposition is an immediate consequence 
of general facts on attractive systems, see e.g. [3] (chs. II and III). 

Proposition 3.1. Let tjo and C,o be two particle configurations such that rjQ < C,o, and let 
P,,o, respectively P^^, be the law of the process starting from rjo, respectively ^q. Then there is 
a coupling Q of P,,o and P^^ (i. e. Q is a measure on the product space, with P^^ as its first 
marginal, and P^^ as the second one) such that 

0--Vt< 6,Vt} = l (3.1) 
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Proof. Being well known that the process corresponding to Lq is attractive, it suffices 
to observe that the flip rates c(x,?7) := D±ri{x) in I± are attractive in the sense that if 
?7(x) = ^(x) = and f] < C, then c{x,ri) < c{x,^), while if t]{x) = ^{x) = 1 and t] < C, then 
c{x,^) < c{x,ri). D 

The analogous nionotonicity property holds for the macroscopic equation. Instead of a direct 
proof we derive the result as a consequence of the nionotonicity of the particle system and 
that it converges to the macroscopic equation. 

Theorem 3.2. Let po, Po be bounded measurable functions from [—1,1] to [0,1] such that 
Poij') < Po(r) for all r G [—1,1], and let p{r,t), respectively p{r,t), be the corresponding 
solution of (2.7) with initial datum Pq , respectively Pq. Then p{r,t) < p{r,t) for all r G [—1, 1] 
and t > 0. 

Proof. Let i^'-^-* and v^^^ be the product probability measures on {0,1}^^ such that 
i'\'^> (^r]{x)) = po{N~^x) and h'^'^\r]{x)) = j5o{N~^x) for all x G Atv- It is well known that a 
coupling A*-^-* of u^^^ and z/*-^-' such that A*^^^{(?7,f/) : r] < f]} = 1 exists. Using Proposition 
3.1 and the notation of Theorem 2.1 we have 

E,w (r/(x, t)) < E^w (r/(x, t)), Vx G A^v, Vt > 0. (3.2) 

From (2.4) we then have that for all t > and for all r G [—1,1], (below [■] denotes the 
integer part) 

p{r,t) = lim E,w(r/([iVr],t)) < lim E^w (^^([iVr], t)) = p{r,t). (3.3) 

u 



4. The macroscopic profile 

We ffist prove that the function p* in the statement of Theorem 2.2 is a stationary solution 
to the Dirichlet problem (1.1) with boundary condition (2.5) or, equivalently, of the integral 
equation (2.7). In fact by requiring that a stationary solution is a linear function we get, 
due to (2.6), that the values of this function at ±1, denoted with u±^ must satisfy 

j(l-n^)=j(l-(l-^-)^)- 



This implies 



2^-^ ~ ^ ■'- K\ „, ^ , -/I „.K\ 



M+ = (l-u_), and =j(l-u+), u+ = -+j(l-^+ 



2 



Solving we get 

in agreement with (2.11). 



J<-^) =3^\ 
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On the other hand, since ■^Pt{r,r') = ^Q^jwPtii^,^') and it satisfies Neumann boundary 
conditions at ±1 we easily see that 

- / Ptir, r'Y'dr' = - {Pt{r, -1) - P^ir, 1)) . 

at j[_i,i] ^ 

Recalhng (from (2.11)) that J = j{l - (p*(l))^') = j(l - (1 - P*(-l))^) we see at once that 
p* satisfies (2.7), which in this case can be written as: 

p*{r) = P,p*{r) + 1(1 - {p*{l)f) J {Ps{r, 1) - Ps{r, -l)}rfs, (4.1) 

for all t > 0. 

We now prove that any solution to the Dirichlet problem converges exponentially fast to p* 
as t — )■ oo. In particular, one has uniqueness of the stationary solution. 

Theorem 4.1. There exist positive constants c, c' so that for any function po G L°°i['-'^i I], [0, 1])| 
the solution p{r, t|po) of the macroscopic equation (2.7) with initial datum p{r, 0) = po{r) sat- 
isfies 

sup |p(r,t|po)-p*(r)|<c'e-'=*. (4.2) 

re[-l,l] 

Proof. Let p{r, t) denote the solution with initial datum p = 1, and p{r^ t) that corresponding 
to initial datum p = 0. From Theorem 3.2 we know that p{r, t) < p{r, t\po) < p{r, t), for any 
initial po. Hence, calling 

w{r,t) := p{r,t) — p{r,t) > 0, w(t) = sup w{r,t) 

~ r6[-l,l] 

it suffices to show that w{t) < c'e~^^ for suitable positive constants c, c' and all t > 0. 

In the proof below c, c, c will denote suitable positive constants (that might depend on the 

model parameter j) whose value may change from line to line. Let 

M±(t):=p(±l,t), M±(t):=p(±l,t), w;±(t):=n±(t)-u±(t)>0. 

/^From (2.7) we see that for all r G [—1, 1], and all t > to > 0, 

w;(r,t) = (P,_,„w;(-,to))(r)-| / f{r,s,t-s)ds, (4.3) 

^ J to 

where 

/(r, s,t-s) := Ps{r, 1) {u+{t - s)^ - u^t - s)^} 

+Ps{r, -1) {(1 - u_{t - s)y' - (1 - u4t - s)f} . (4.4) 

Interchanging particles and holes, one can couple at once the evolutions starting from the 
configurations f] = I {all occupied sites) and rj = {all empty sites) so that f]{x,t) = 
1 — r]{—x,t). Therefore, by the same argument as in the proof of Theorem 3.2 one has 
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p{r,t) = 1 — p{—r,t) for all r and all t. In particular w{—r,t) = w{r,t), u±(t) = 1 — u^it) 
and w+{t) = w^(t) for all t. (Still from Theorem 3.2 we see that w{r,t) and so also w(t) 
decrease in t.) Of course w{r, 0) = 1 for all r. 
In particular, we may rewrite (4.3) with to = as 

w(r, t) = 1 - ^ /" [P,{r, 1) + Psir, -1)]«;(1, t - s)h{t - s)ds (4.5) 

where 

K-l 

h{t -s):=J2 ^+it - s)''~'~'u+it - sY (4.6) 

and where we have used that for any integer K > 1, 

a'^ - h'^ = [a-h)Y^ h'a'''^~', a > 6 > 0. (4.7) 

1=0 

Also, from (4.6) and the monotonicity properties we see that 



h ■= p*{ir~' < h{t) <b + K-l=:cK. (4.8) 

imes. To this end we introduce the kernel operat( 

Ki'^fir) = - [ Ps{ry)f{r')dr\ f G C([-l, 1], 

following integral equation: 



The proof will use local times. To this end we introduce the kernel operators Kg , e > 0: 



'[-l,-l+e]U[l-e,l] 

In particular Kg f{r) = Kg f{—r) for all r G [—1,1]. Let w^'^^ be the solution to the 



w^'\r,t) = 1-- [ {Ki'^w^'\-,t - s)){r)h{t - s)ds. 
2 Jo 

We shall next prove that for all T > 0, 

lim sup sup \w{r,t) — w^^\r,t)\ = 0. (4.9) 

<^^o re[-l,l] 0<t<T 

Calling 

ilj{r,t) = \w{r,t)-w^'\r,t)\, ^(t) = sup tlj{r,t) (4.10) 

and using (4.8), we can write 

/ {iKi'^w^'\;t-s))ir)-{Ps{r,l) + Ps{r,-l)}wil,t-s)}h{t-s)ds < ce 
Jo 

+CK f <-[ \Ps{r,y)-Ps{r,l) + Ps{-r,y)-Ps{-r,l)\dy\\w^'\y,t-s)\ds 
+CK J {Psir, 1) + Psi-r, l)}*(t - s)ds. (4.11) 
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Using that for all y G [1 — e, 1], r G [—1, 1] 

\Ps{r,y) - Ps{r, 1)| < c^, Ws G [e,t] (4.12) 



we see that the second term on the r.h.s. of (4.11) is bounded above by 

for suitable constants c, c'. We then easily get 

^{r,t) < civ^+ C2 / ^(s)cis (4.13) 

Jo 

for suitable constants ci,C2. By the Gronwall inequality we conclude (4.9). 

We now estimate w^^\ Let {Bt} be a standard Brownian motion with reflecting b.c. at 
±1, with Fr denoting its law when Bq = r (and corresponding expectations denoted by E^.). 
Then 

where 



ip,{B,t-s)=<P,iB)h{t-s), 0,(r) = ^l[i_,,i](|r|), rG [-1,1]. (4.15) 



By (4.8) 



w^^'ir 



,t)<Er(e-''^oMBs)dsY ^^-^Q^ 



For < t < t we write 



We shall prove below that taking i sufficiently small, we can take a < 1 so that for all e > 

sup EJe-^^o^^(^=)'^A < 1 - a (4.18) 

r-e[-l,l] ^ ^ 

From (4.18) and (4.17) we then get 

\w^'\r,t)\ < {1 - aY'^^ (4.19) 

([a] the integer part of a) which then concludes the proof of the theorem. 

Proof of (4.18). Let T = inf{t > : \Bt\ = 1}. We then have 

Er(e-''^oMBs)ds\ < Er(l{T<i/2}e~^^^^'^^'^'^') +f'r{T >i/2) (4.20) 
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and write 

< P,(T<t72)Ei(e-''^o^''^^(^^)'^')) (4.21) 



where we also used that Ei(e-^i'o^''^=(^^)'^^) = E_i(e~''^o^''^^(-B-)'^'^) by symmetry. 
By Taylor expansion 



nt/2 

El (e-''^o^' MBs)ds^ < 1 - 6Ei ( / (pe{Bs)dsj + 6 (4.22) 

where 

6 = (|^)' / dh f' dh [ Pt,{l,yi)Pt,{yuy2)dyidy2 (4.23) 

^e Jo ^0 ^|»/i|G[l-e,l],|j/2|e[l-e,l] 

But, from (2.8)-(2.9) we see that 

sup Ps{x,y)<c— (4.24) 

x,yg[-l,l] yS 

so that for i small enough we get 

6 < ci/2 (4.25) 

for suitable constant c. Using again (2.8)-(2.9), we see at once that a positive constant c 
can be taken so that for all i small, and all e > 

bEi( (j),{Bs)ds) > c^/t/2. (4.26) 

From (4.22), (4.25) and (4.26) we then get for t small (with possibly different constant c). 

El (e-^^o'" MBs)ds^ < I _ c^/t/2. (4.27) 

By (4.20) and (4.21) we then have 

Er(e-''^o'^'(B^)'iA <[1- c^/t/2]¥r{T < i/2) + P^(T > t/2) < 1 - a (4.28) 



with 



a = inf P^(r < t/2)c^/tj2. (4.29) 

re[-l,l] 
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5. Proof of Theorem 2.2 

The proof is a direct consequence of the following three facts, (i) For any t > and any 
integer k > 1 



lim max max 



Ar^oo^g{0,l}A]v (xi,..,Xk)eA''-^ 



E^{Y[vix^,t))-Y[E^{r]{x,,t)) =0 (5.1 



(ii) For any t > 



lim max max 



E,{r]{x,t))-p{x/N,t\r]) =0 (5.2) 



m 



lim sup ||p(-,t|po)-p*(-)lloo = (5.3) 

*^~PoeL°°{[-l,l];[0,l]) 



(i) and (ii) are proved in [2]-[l], (5.3) is proved in Theorem 4.1. 
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